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ABSTRACT 

CFHT Megacam data in (g',r',i') are used to obtain deep, wide-field photometry of the globular 
cluster system (CCS) around M87. A total of 6200 GCs brighter than i' = 23.0 (roughly equivalent 
to Mi — —8.5) are included in the study, essentially containing almost the entire bright half of 
the total GC population in the galaxy. The classic bimodal metal-poor and metal-rich sequences of 
GCs show up clearly; while the spatial distribution of the GCs can be traced detectably outward to 
R gc ~ 100 kpc and perhaps further, the blue, metal-poor subpopulation is very much more spatially 
extended than the red subpopulation. Both the red and blue GC subsystems have radial metallicity 
gradients, where mean heavy-element abundance scales with projected galactocentric distance as 
Z ~ R~ 012 (blue) and i?~ 017 (red). The blue sequence exhibits a strongly significant mass/metallicity 
relation (MMR) in which the mean metallicity gradually increases with cluster luminosity as Z ~ 
£0.25±o.05 f or j-jjg l um i n osity range Mi < —10 and the assumption of a constant M/L. However, this 
relation is also clearly nonlinear: fainter than this level, the sequence is more nearly vertical. This 
mass/metallicity trend can be understood as the result of self-enrichment within the most massive 
metal-poor GCs during their formation. The red sequence formally exhibits a negatively sloped MMR, 
but the numerical solutions and tests show that this red-GC slope is not very significant. In giant 
elliptical galaxies, the red GCs are likely to represent a broad composite population formed during 
several major starbursts. If so, the red sequence might display a population-based MMR that could 
in principle be either positive or negative. 
Subject headings: galaxies: elliptical and lenticular, cD 



1. INTRODUCTION 

By far the largest collection of globular clusters (GCs) 
within any one galaxy in the local universe is in M87, the 
cD giant in the Virgo cluster. Its total GC population 
is N ~ 14000 (|McLaughlin et alJ H99l iTamura et all 
120061 ). roughly a third of the GCS population in all the 
Virgo gal axies combined (s ee the Virgo Cluster Survey 
results of iPeng et alj|200a ). In addition, because both 
the foreground reddening and field contamination are 
low, the M87 GCS can be traced outward to extremely 
large projected galac tocentric distance, 100 kpc or more 
(jTamura et al.l 12006) . This combination of characteris- 
tics makes M87 a unique testbed for systematics of GC 
systems. 

A particular topic of recent interest is the correlation 
between GC metallicity and luminosity (or mass), a trend 
which has been found to show up on the blue, metal- 
poor sequence (lHarris et al.l 120061: IStrader et aLl | 2006c 
Mieske et al J I2006t ISpitler et alJl2006t|Forte et aljl2007t 
Wehner et alJl2008tlBassino et alJl200alHarrisl2009l and 
references cited there). In brief, the mean GC metallic- 
ity along the blue sequence is seen to increase steadily, 
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though perhaps not linearly, with increasing cluster 
mass. So far, the most likely physical explanation that 
has been advanced for the origin of this effect is some 
form of self-enrichment during the protocluster epoch 
(lHarris et alj 120061: IMieske et alJfeOOfl : IStrader fc Smith! 
120081: iBailin & Harrisl2009| r The essential idea of this in- 
terpretation is that any self-enrichment should be more 
efficient within the more massive, deeper potential wells 
of the larger clusters. 

M87 has turned into something of a flashpoint for 
discussion of this mass/metallicity relation (MMR). 
IStrader et al.l (|2006D first claimed that the mean GC 
color along the entire blue sequence increased roughly 
linearly with absolute magnitude, which would corre- 
spond to a simple power-law scaling of heavy-element 
abundance with cluster mass, Z ~ M p , for a con- 
stant mass-to-light ratio. Their measurements con- 
sisted of simple fixed-aperture photometry of ~ 1000 
GCs around M87 from the single-orbit exposures in 
(</, z') of the galaxy that were taken as part of the 
Virgo Cluster Survey (|C6te et al.l I2004D . IMieske et all 
(2006) used the same material in the full Virgo Clus- 
ter Survey analysis in which the cluster photometry was 
done with PSF-convolved King profiles. They found 
much the same trend though it was less clear that the 
MMR was as stra ightforward as a li near fit in the color- 
magnitude plane. iForte et al.l (|2007l ) used ground-based 
C, T\ photometry with standard PSF-fitting photometry, 
and found a roughly linear correlation of ( C — T\) color 
versus magnitude along the blue sequence. I Waters et al.l 
(2009), on other hand, find no net MMR along the blue 
sequence (that is, they find a vertical sequence of magni- 
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tude versus color) by using size-corrected aperture pho- 
tometry of about 2000 clusters from deeper HST/ACS 
image s of M87 in (V,I). Most recently, iPeng et all 
(2009) have analyzed the same ultra-deep HST/ACS 
data with the PSF-conyolved King profile technique de- 
scribed in Uordan et al.l (|2009f ) to find that a significant 
color/magnitude slope exists for the upper half of the GC 
magnitude range (roughly, for the clusters brighter than 
the turnover point of their luminosity function). Thus 
five different measurement procedures from four partially 
different datasets have yielded results that are largely in 
mutual agreement, but also show points of disagreement 
that are still of some importance and need to be settled. 

For comparison, the largest GC database used specif- 
ically to analyze the MMR is the sample of ~ 10, 000 
clust ers in six giant Brightest Cluster Galaxies (jHarrisI 
2009). There, the photometry of the individual GCs is 
obtained from PSF-convolved King-model profiles and is 
taken in the metallicity-sensitive (B — I) index and with 
deep exposures. 1 For clusters brighter than the GCLF 
turnover the data have high S/N, and the total sample 
is so large that the GCLF can be traced upward to the 
most luminous known GCs at L > 10 7 Lq. Extensive ev- 
idence is presented that the MMR has a nonlinear form: 
the blue GC sequence is nearly vertical (i.e., no trend of 
GC color versus luminosity) for L < 5 x 10 5 L Q and then 
curves toward redder color at higher luminosity. 

The quant ita tive self-enrichment model of 
iBailin fc Harris! (|2009f ) reproduces this nonlinear 
MMR form fairly accurately. The key physical point 
in the model is that protoclusters less massive than 
about 1O 6 M do not have deep enough potential wells 
to hold back much of the first round of SNII-enriched 
gas, so for these lower-mass GCs there should be no 
MMR on either the blue or red sequence. Conversely, 
for the higher-mass GCs with their deep potential wells, 
self-enrichment is almost automatically required. A 
unique prediction of this particular model is that both 
the red and blue GC sequences should exhibit a MMR 
at the top end, though it should be less noticeable for 
the red sequence because it starts from a much higher 
level of pre-enrichment. 

Settling the issue on observational grounds of the de- 
tailed form of the MMR (and indeed whether or not it ex- 
ists in all galaxies) is of some urgency because confusion 
over the form of the correlation has already arisen in the 
literature. Yet it is likely to bear quite directly on under- 
standing the formation of massive GCs and their early 
star f ormation history (iHarris et al.l 2006; Mieske et al.l 
[20061 IStrader fc Smith! [20081 : IBailin fc Harridl2009n~ 

In this paper I discuss new photometry of the M87 
GCs obtained with CFHT/Megacam. The data com- 

1 In this and most previous MMR work, the fundamental but 
well justified assumption is made that the spread of broadband 
optical color indices directly represents differences in metallicity, 
with [Fe/H] increasing monotonically as the color index gets red- 
der. The effects of other parameters on their broadband colors 
(such as cluster age or a— abundance) generate far more subtle 
changes in color. Convincing measurements of age or a differ- 
ences between old GCs traditionally require careful analysis of 
high— S/N integrated spe ctra and are challen g ing to uncover eve n 
then; see, fo r example. IPuzia e t al. (200^); IPene et al.l (2006); 
IStrader et all pOOfih : IPessev et alj H2008I) among others, as well 
as compari sons of contempora ry simple stellar population (SSP) 
models (e.g lPessev et al.|12008t) . 



prise a large enough GC sample (N ~ 10 4 ) to analyze 
the MMR to deep limits in addition to correlation with 
galactocentric distance. The sample size provides the 
obvious advantage that features such as the MMR, the 
spatial distribution, or the GCLF can be studied within 
a single galaxy without having to put together composite 
samples. The results discussed in the following sections 
show well defined blue and red GC sequences extending 
to high mass, and exhibit a nonlinear MMR quite simi- 
lar to the ones in the six-galaxy BCG sample. Though 
M87 alone can hardly settle the question of the univer- 
sality of the MMR, these results reinforce the need to 
investigate this intriguing effect further. In the present 
paper, I concentrate on the various features of the color 
and metallicity distributions. 

In this discussion I adopt a distance modulus (m — 
M)n = 31.0 for M87 and the core of the Virgo cluster 
(jJerien et al.l 120041 : ICaldwelll 120061 : 1 Williams et al.ll2007h 
and a foreground absorption and reddening Aj = 0.04, 
Ev-i = 0.03. In the SDSS filter system the correspond- 
ing values relevant here are Ai = 0.04, E g -i = 0.04, 
E g ^ r = 0.02. 

2. THE DATA SAMPLE 

The 36 CCDs making up the CFHT Megacam array 
provide a 1° square field of view and the camera is nor- 
mally used with the u'g'r'i' z' filter set (in what follows I 
will drop the accents and refer to the magnitudes simply 
as ugriz). The image scale is 0.187"/px and the entire 
field makes up a 20K x 20K array. The data used in this 
paper consist of a single set of homogeneous exposures 
centered on M87 taken in gri during 2004 March and 
April (Harris, PI), as listed in Table 1. While this ma- 
terial will eventually be surpassed by the much more ex- 
tensive CFHT/Megacam Next- Generation Virgo Cluster 
Survey by Ferrarese et al. (see www.astrosci.ca/NGVS), 
much can already be extracted from the earlier set of 
M87-centered exposures. 

Individual exposures in each filter were taken with a 
"small dither" pattern sufficient to splice over the small 
gaps between the individual CCDs in the array, but the 
much larger gaps between the top two ranks of CCDs 
and the bottom two ranks are still present (standard 
large-dither patterns are available for the camera that 
would cross these gaps, but these make the astromet- 
ric re-registration more difficult). These remaining gaps 
show up as two ~ 350— pixel- wide blank E/W stripes 
in the array, making up about 3.5% of its total area on 
the sky. The hour-long total exposure times in each fil- 
ter were chosen to reach closely similar limiting absolute 
magnitudes for the globular clusters. The combined im- 
ages in each filter are the ones produced through the 
Megapipe facility of CADC, and were downloaded from 
the CADC Archive ready for photometry. 

The GCs around M87 are seen against a population of 
field objects including foreground Milky Way stars and 
background galaxies. Because the Virgo field is at high 
latitude, the contamination from stars is minimal and the 
main contaminants are small, faint galaxies. However, 
much of the field galaxy population can be removed by 
morphology because the GCs themselves appear starlike 
on these ground-based images. At the 16-Mpc dist ance of 
Virgo , an average GC half-light diameter of 5 pc {Harris 
1996) is equivalent to 0.06", which in turn corresponds 
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to a FWHM of 0.04" for an average Ki ng-model concen- 
tration log (r t /r c ) = 1.5 (|Larseniri999f) . This scale size 
is thus 20 times smaller than the 0.8" FWHM seeing on 
the combined Megacam frames in this dataset, making 
any such objects indistinguishable from stars. Exten- 
sive quantitative tests of PSF-convolved globular cluster 
profiles for distant galaxies such as these show explicitly 
that the intrinsic profile width falls below the threshold 
of measurement if i ts FWHM is l ess than 10 percent of 
the PSF width (e.g. ILarsedfl999l IHarrisll20l)9h . In prac- 
tice, what this means for the present data is that the 
GCs can be accurately measured either through conven- 
tional fixed-aperture or PSF-fitting photometry. Because 
the crowding levels are minimal, aperture photometry is 
adopted here (see below), with an aperture radius set 
at its optimal value similar to the PSF FWHM. Use of 
aperture photometry also alleviates concerns over minor 
changes in the PSF shape across the field. 

As another basis of comparison, it is worth noting that 
imaging of the M87 clusters with ground-based image 
quality of 0.8" is almost exactly analogous to imaging 
clusters in the Coma galaxies at d = 100 Mpc with HST 
image quality of 0.1". Explicit test s of PSF-convolved 
GC profiles for the Coma galaxies (Harris 
with deep HST data verify that the intrinsic sizes of the 
GCs are too small to affect their integrated magnitudes 
from either PSF-fitting or fixed-aperture photometry at 
levels larger than > 0.02 mag. Even further, when these 
second-order magnitude corrections are differenced be- 
tween two filters (g — i or g — r in this case), the effects of 
aperture size corrections on color are even smaller third- 
order ones that can be completely neglected. Similar 
comments apply to the (C — Ij ) ground-based photome- 
try of M87 and NGC 1399 bv lForte et all (pOOl in which 
the seeing was 1.0" — 1.6". 

The first step in the photometry was to find candi- 
date objects on all three master images and to eliminate 
the majority of nonstella r objects. The SExtractor code 
(|Bertin fc Arnoutslll996f) was run to do this, and plots 
were generated of the SE effective radius and FWHM ver- 
sus aperture magnitude, as shown in Figure □ Objects 
falling outside the generous "stellar" boxes marked in the 
figures are either very bright, saturated stars; obviously 
nonstellar extended galaxies; or extremely faint objects, 
all of which could be eliminated from the detection lists. 

Next, the SE-cullcd finding lists defined separately in 
g, r, i were compared, and any objects found on all three 
filters and coordinate-matched to within ±1.5 px were 
kept for the final photometric measurement run. These 
matched objects were measured through daophot with 
an aperture radius of 3.5 px (0.65", equivalent to a linear 
size of 50 pc at the distance of M87). Crowding in this 
high-latitude field is negligible, and comparisons of the 
fixed-aperture magnitudes with PSF-fitted magnitudes 
showed no differences in the resulting color-magnitude 
diagrams. 

3. ASTROMETRIC AND PHOTOMETRIC 
CALIBRATION 

The Megacam stacked images in the CADC Archive 
have internally homogeneous astrometry that is already 
close to the most commonly used systems (such as SDSS, 
2MASS, USNO), but to ensure that the data were on 
a specific astrometric system as closely as possible, the 



USNO-B1 system was selected. The USNO-B1 cata- 
log was searched for all objects within the Megacam 
field, and the positions of these objects were matched up 
with the (a, 6) values returned directly from the Mega- 
cam image headers. This matchup yielded 932 stars in 
common to within 1". Small zero-point offsets equal to 
Aa = (0.053±0.008)", AS = (0.254±0.008)" in the sense 
(USNO-CFHT) were found, but no significant higher- 
order terms. Figure [2] shows the residuals after the na- 
tive Megacam astrometric coordinates were shifted onto 
the USNO-B1 system. The object-to-object rms scat- 
ter is ±0.23" across the entire field, with no systematic 
deviations larger than ±0.05". 

The Megacam ugriz magnitudes are known to differ 
slightly from the Sloan Digital Sky Survey scales (see 
the Megacam webpages for updated calibration informa- 
tion). However, for the purposes of more direct compar- 
ison with other GCS data in the literature, it is useful 
to put the present M87 data onto the SDSS scale. The 
SDSS DR5 catalog was searched for objects in the M87 
field and compared with the Megacam object lists, with 
the result that matches were found for about 1000 stars 
brighter than i = 21.0. These were used directly to set 
the photometric zeropoints for gri in the M87 field. The 
mean differences Am = m(CFHT) — m(SDSS) were 
found to be A{g) = 0.254 ±0.004, A(r) = 0.407 ±0.008, 
and A(i) = 0.412 ± 0.006, where the CFHT magnitudes 
are the ones defined simply using the zeropoints in the 
image headers. 

A secondary advantage of tying these data closely to a 
widely used standard system is that the intrinsic colors 
of globular clusters in the (g — r) and (g — i) indices have 
not previously been well established. The Virgo Clus- 
ter Su rvey data are in (q — z) through th e HST/ ACS 



filters ([Peng et all 120061 : iMieske et al.ll2006h . Previously 
publish ed studies of GCSs f o r various galax i es in ( q — i) 
include iForbes et all (12004D . iWehner et all (|2008f) . and 
iCockcroft et alj (|2009f ) all taken with Gemini/GMOS, 
but in all of these the absolute photometric calibrations 
are preliminary and not likely to be accurate to better 
than ±0.1 mag (see these papers for detailed discussion). 
M87 offers an excellent opportunity to set up a fiducial 
color scale for GCs, not just because a direct matchup to 
the SDSS catalog is possible, but also because its fore- 
ground reddening is almost negligible and the statistical 
size of the GC population allows the mean colors of the 
red and blue sequences to be defined internally to ±0.02 
mag or better (see below). 



4. DEFINITION OF THE GLOBULAR CLUSTER 
SYSTEM 

In Figure [3] the final data for 20189 starlike or near- 
starlike objects are shown, plotted in (i, g — r) and (i, g — 
i). The detection rate for objects within the color range 
(g — r) < 1.0 is highly complete for i < 23.5; in the 
following analysis, the main conclusions rely even more 
conservatively on the interval i < 23.0 and are unaffected 
by any incompleteness. 

Several clearly distinct components show up. The GCs 
that we are after lie on the two nearly vertical, narrow 
sequences at [g — i) ~ 0.8, 1.05 while the populations 
bracketing these on the blue and red sides are the ex- 
pected loci of foreground Milky Way halo and disk stars. 
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Any objects falling in the regions 

{g - r) < 0.43, (g - r) > 0.95 
(g - i) < 0.59, (g - i) > 1.40 
i < 18.5 

(the boundaries of which are marked by the boxes in 
Fig. [3]) are thus rejected as field contamination. This 
final culling leaves 11946 objects that make up our final 
"best" GC candidate list. 

The distribution of the GC candidate colors by spatial 
location is shown in Figure |U Here, the projected dis- 
tance R from the center of M87 is expressed in units of 
Reff — 1-58', the effective radius of the galaxy's spheroid 
light. A well known though still striking feature of the 
M87 GC system is its huge radial extent: the clusters 
extend detectably beyond even R ~ 150 kpc or 20 R e ff 
(see next section). The first hints of its huge GCS halo 
were already ev i dent in the wide-field starcount work of 
IHarris fc Smith! (fl976h and m ore clearly outlined with 
the CCD-based ph otometry of IMcLaughlin et alj (|1994l ) 
and lTamura et all (120061). For comparison , howev er, the 
V— magnitude data of IMcLaughlin et al.1 (|1994D could 
be used to trace the GCS out only to the radial lim- 
its of their single-CC D field at R(max) ~ 8' ~ 40 kpc. 
iTamura et al.l (|2006l ) used a combination of four Subaru 
Suprime-Cam pointings in three colors (BVI) to make 
up a surveyed region of 136' x 27', about the same total 
area as the single Megacam field but oriented along the 
E/W direction. M87 itself is in the SW corner of their 
westernmost SCam pointing, while another Virgo giant 
NGC 4552 is in their third field pointing. 

An equally striking feature of the spatial distribution 
from Fig. [4] is the difference between the blue and red 
GC subpopulations. The red sequence virtually disap- 
pears into the general spread of field contamination for 
R ^ 6 ~ 8i2 e // (~ 12' or 50 kpc) while the blue sequence 
continues detectably outward much further. Two other 
illustrations of this are shown in Figures [5] and O As 
a first approximation we define the blue GCs as those 
in the color interval 0.59 < (g — i) < 0.93, and the red 
GCs as those in 0.93 < (g - i) < 1.40. In Fig. El the 
spatia l pos itions of each type are shown. ITamura et alJ 
(|2006h and lForte et~aT] (|2007h find exactly the same phe- 
nomenon and these papers provide a more detailed dis- 
cussion of the spatial distributions of each component. 
In Fig. [6l the CMDs in (i, g — i) are displayed for five 
radial zones. Whereas the blue and red sequences are 
almost equally populated for the innermost zone R < 3', 
the red sequence dies away more rapidly outward, leaving 
the blue sequence to dominate at large radii. 

5. RADIAL STRUCTURE 

The M87 GCS is so populous and extended that the 
level of background contamination (that is, the field con- 
tamination still present in the list of starlike objects, as- 
sumed constant across the field) is hard to determine 
precisely. The contamination is completely unimpor- 
tant at small radii, but at the largest radii of the sur- 
vey it determines the effective limit beyond which the 
GCS can be traced. We have no truly remote "control 
field" adjacent to M87 for measuring <7b, and ultimately 
a definitive measurement of the background field level 



and the outermost radial profile of the M87 GCS will 
have to await the still wider-field work from the Next- 
Generation Virgo Survey. For the present, the outer 
parts of the single Megacam field itself must be used to 
gauge <j\>. A reasonable estimate can be made by noting 
that for R > l2R e ff, the areal density of red objects 
(i < 23.0 and 0.93 < (g — i) < 1.4) is nearly uniform 
at (7 = 0.36 arcmin -2 , accounting for the edges of the 
frame and the two blank stripes in the CCD array. The 
color range enclosing the blue GC sequence is narrower 
(0.59 < (g — i) < 0.93), so if the color distribution of the 
contaminants is roughly uniform over this range, the level 
of the blue background will be proportionately smaller, 
&bkdg(blue) ~ 0.26 arcmin -2 , for the same magnitude 
range i < 23.0. (See also Harris 2009 for more extensive 
tests and discussion of the background color distribution 
in similar high- latitude fields.) Thus for the present pur- 
poses I adopt a total (Jbkgd — 0.6 arcmin -2 , with a true 
external uncertainty likely to be near ±0.1 arcmin -2 . 

In the final list of GC candidates defined above, there 
are 8313 objects brighter than i = 23.0. Subtracting 
off the background, this leaves a total GC population 
N GC ~ 6200 ± 350 brighter than Mj ~ -8.5, which is 
about 0.2 mag brighter than the turnover (peak) point 
of the luminosity function. 

The actual radial profile of both the blue and red GC 
subsystems can be quantified by tracing the projected 
surface density a c i directly as a function of radius. The 
observed number density of all the candidate GCs is 
by definition a = a c i + Ub- The resulting background- 
subtracted data for i < 23.0 are listed in Table O and 
shown in Figure [7] The columns in the Table give the in- 
ner and outer radii of each annulus, the mean radius, the 
calculated total number of clusters in the annular zone 
after background subtraction, and the density a c i. The 
conventional log/log form (log o c i vs. log R) is shown 
in the left panel, indicating that both the blue and red 
components fall off progressively more steeply with in- 
creasing radius. A standard Hubble profile, which is a 
power law with a core radius Rq, 

a d (R) = + (l) 

provides a reasonable fit to either component. Fits of 
this model to the data (shown in Fig. [7|) for the blue 
GCs give (cto = 66 arcmin -2 , Rq = 2.0', a = —1.8); 
while for the red GCs, (oo = 150 arcmin -2 , i?o = 1-2', 
a = —2.1). These are equivalent to simple power laws 
a ~ i? - ™ where the exponent n steepens continuously 
with radius, 

n ( R ) = T\~TJL\ ' (2) 

Both components also follow a deVaucouleurs-law form 
(log a c i vs. i? 1 / 4 ) fairly well over their entire observed 
run. In this form the best-fit equations describing each 
one are 

loga c i(blue) = (3.154 ± 0.067) - (1.760 ± 0.045)(i?/i? e// ) 1/ t3.) 
\oga c i(red) = (3.874 ± 0.054) - (2.439 ± 0.037) {R/ R ef f) 1 ^^) 

where again the a's are in units of arcmin -2 and the 
slopes are normalized to the Reff of the halo light. 
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(These lines should, however, not be extrapolated in- 
ward where the density di stribution flattens of f more 
strongly.) For comparison, iTamura et all (|2006f) found 
slopes (-1.47±0.03)(blue) and (-2.44 ± 0.06)(red). As 
they discuss more extensively, the integrated light of 
the M87 halo has a fl^-law slope of -2.36 ± 0.09, 
strongly suggesting that the more metal-rich GC pop- 
ulation can be directly linked to the similarly metal-rich 
population of the field halo stars. Much other discus- 
sion of similar results from several other galaxies re- 
inforcing this interpretation and conn ecting with their 
evoluti onary histories c an be fou nd inlLee et alj (119981 . 
20081): [Harris fc Harris! (I2002T): IBassino et all 12006) ; 
Forte et all (|2007l) : lShin fc Kawatal (|2009ft among others, 
and will not be repeated here. 

The behavior of the deVaucouleurs-type profile fit it- 
self, near its outer edge, provides a useful consistency 
check of our adopted backgrounds cr& (blue, red). If we 
have chosen the background nearly correctly, then the 
resulting er c ; profile can plausibly be claimed (as is the 
case for most other elliptical galaxies) to follow a nearly 
straight line in the log <j c i vs. r 1 / 4 plane all the way 
to the outer edge of the data. An overestimate of o\, 
would cause the background-subtracted profile to drop 
suddenly and unrealistically to zero at the edge, whereas 
an underestimate would cause a c \ to flatten off artificially 
at large radius. Neither one is the case (Figure 7b), so the 
choice of background levels appears to be self-consistent 
within their stated internal uncertainties. 

The a c i data indicate that the blue GCs exist in de- 
tectable numbers out to at least R = 28' or ~ 130 kpc. 
A final c heck of the present re sults comes from compari- 
son with lTamura et aLl (|2006h . whose SCam data extend 
along a line towards NGC 4552 out to 120' (more than 
500 kpc) from M87. They find the blue-sequence GCs 
to extend detectably to R = 30' (see their Figures 5 and 
7), and present a model fit to the combined M87/N4552 
counts that suggests the M87 GCs could even extend to 
100' and more. Clearly, we have not yet traced the full 
extent of this system. 

6. METALLICITY GRADIENTS 

A second-order feature of the radial color distributions 
(Fig. |4j) is a progressive change in the mean color of both 
the red and blue GCs with radius: in other words, both 
components show a metallicity gradient. This feature is 
clearest for R < 8R e ff or about 60 kpc; beyond that 
point, the red clusters are almost absent and the blue 
clusters show no significant change in mean color. For 
the region inward of 8i? e //, power-law fits yield 

{g - i)(blue) = 0.813 - 0.0251og(i?/i? e// ) , (5) 
(g-i)(red) = 1.095 - 0.0351og(i?/i? e// ) . (6) 

where both slopes are uncertain by 10 percen t . Th e 
single-stellar-population models of iMarastonl (2005), 
which provide integrated color indices for a variety of 
broadband colors including the SDSS filters, show that 
A(g-i)/A[Z/H] = 0.21±0.05 over the metallicity range 
[Z/H] < —0.2 that covers most GCs. The color slopes 
then convert into metallicity gradients Z ~ R 1 where 
7 = 0.12 ± 0.02 (blue) and 0.17 ± 0.03 (red). 

Comparable data for the GCSs i n several othe r large 
galaxies are now in the literature. iGeisler et al.l (|1996f ) 



found 7 = 0.15±0.03 (blue), 0. 12±0.06 (red) in th e other 
Virgo supergiant NGC 4472. iForte et all (j2001h found 
7 ~ 0.2 for the blue GCs in two Fornax ellipt icals, NGC 
1399 and 1427. Most recently, lHarrisI (j2009l ) derived a 
mean 7 = 0.10 ± 0.02 for both the blue and red GCs in 
six giant Brightest Cluster Galaxies. All these results are 
similar, and begin to present a very consistent picture: 
the mean GC metallicity in these large galaxies depends 
on spatial location, becoming more enriched the further 
inward it resides in the halo. 

The widespread presence of metallicity gradients, al- 
beit shallow, indicates that both the major epochs of GC 
formation, metal-poor and metal-rich, had a significant 
component of dissipative collapse of protogalactic gas 
and in situ formation which led to higher enrichment lev- 
els deeper in the potential wells of these large galaxies. 
Although gradients in the range 7 ~ 0.1—0.4 for the dom- 
inant metal-rich field-star component have been found in 
many ellipticals within their central bulges (R < i? e //), 
the globular cluster data indicate that dissipative for- 
mation extended much further outward into their halos, 
to ~ 5R e ff and beyond. What is perhaps equally im- 
portant for understanding galaxy formation is that the 
metal-poor component, which is all but undetectable in 
the integrated field-star light, also has an inbuilt metal- 
licity gradient to the same degree. This result makes 
it unlikely that the metal-poor stars all formed in pre- 
galactic dwarfs that only later began to merge into big- 
ger galaxies. Instead, the first major metal-poor round 
of star formation must have taken place while the ini- 
tial stages of hier archical merging were already underway 
(see [Harris 200j| for further discussion and references) 

7. BIMODAL SEQUENCE FITTING 

We now look more closely at the blue and red GC 
sequences in the color-magnitude diagrams. The goal 
is to characterize each of the two sequences through 
their mean colors as a function of magnitude (the 
MMR) and their internal color (metallicity) dispersions. 
To define the MMR (or lack of one) we adopt the 
sam e basic approach as used in several previous studies 
Ce.g.lHarris et al 1120061 iMieske et al.li 2006: We hner et all 
120081: IBassino et al.l 120081: Harris! 120091 among others): 
the CMD is divided into narrow magnitude bins, and 
within each bin the color distribution of the GCs is fitted 
by a bimodal Gaussian model. Because each magnitude 
bin is treated independently, no a priori assumptions are 
made about the particular form of the MMR, whether 
linear or otherwise. 

For the multimodal fitting, the code RMIX is used here 
(for exten sive descriptions of RMIX applie d to this situ- 
ation, see iWehner et al.ll2008t lHamsH2009f ) . The default 
approach - facilitated by the sheer statistical size of the 
M87 GC population - is to allow the code to solve freely 
for all five basic parameters in a Gaussian bimodal dis- 
tribution: the mean colors of the two modes fi2 (blue, 
red), their dispersions cti, 02 (blue, red), and the relative 
numbers (proportions) pxjpi of objects in each mode. 
The MMR along either of the sequences is then defined 
by any trend in the mean fj, with magnitude. 

The bimodal fits were done for objects within the zone 
R < 6.5i? e //, i.e. within the part of the M87 halo where 
both blue and red sequences are definitely present and in 
roughly equal numbers (Fig. [4]) . We restrict the sample 
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to objects in the color range 0.65 < {g — i) < 1-30, and as 
a preliminary step, the shallow metallicity gradient de- 
scribed in the previous section is removed by normalizing 
the color indices to the values that all the clusters would 
have if they were at R = R e ff- That is, for the blue- 
sequence GCs (defined as those objects with (g — i) < 
0.95), we define (g-i) norm = (g-i)+0.0245\og(R/R eff ), 
and for the red-sequence GCs ((g — i) > 0.95), we define 

{g ~ i)norm = (.9 - *) + 0.03451og(i?/i? e// ) . 

The results of these bimodal fits in both CMDs are 
summarized in Table 2, 3, and 4. Four separate runs were 
done with the following conditions, which were chosen to 
allow additional tests of the reality of any deduced MMR 
along both of the sequences: 

• The data were divided into 0.25-magnitude bins 
in i, and within each bin all five free parameters 
(Mi; A*2> cr l) cr 2 ) Pi/p2) defining the bimodal Gaus- 
sian fit were freely determined by the fit (i.e., un- 
constrained). The exceptions to this were at the 
bright end, where for i < 20.5 half-magnitude bins 
were used to capture enough objects for a statisti- 
cally significant fit. In addition, preliminary runs 
showed that for the four brightest bins (i < 21.0) 
the fully unconstrained fits tended to be unsta- 
ble, so for these, the intrinsic widths of the two 
sequences were constrained to be equal to their av- 
erage values on the lower bins, ci = 0.060 and 
(72 = 0.10. The results for this run are summa- 
rized in Table 2. The columns list the i— magnitude 
range defining each bin; the number of objects 
in the bin; the best-fit blue-sequence parameters 
(Hi,cri) and their uncertainties; the red-sequence 
parameters (^2,02) and their uncertainties; and 
the best-fit proportion p\ (the fraction of the to- 
tal population belonging to the blue sequence) . By 
definition P2 = 1 — pi ■ 

• The same conditions as above were used, except 
that the red sequence was constrained to have 
a mean color [12 = 1.070, the average obtained 
from the lower-luminosity bins in the first (uncon- 
strained) run. The color of the blue sequence was 
determined freely from the fit. The results for this 
run are summarized in Table 3. 

• The same conditions were used, except that next 
the blue sequence was pinned at a mean color 
fix = 0.80, the average obtained from the lower- 
luminosity bins in the first (unconstrained) run. 
The color of the red sequence was determined freely 
from the fit. The results for this run are summa- 
rized in Table 4. 

• The mean colors of both sequences were pinned 
(^1 = 0.80, \i2 = 1-07) and only the ct's and pro- 
portions were solved for. 

As an illustration of the overall quality of the funda- 
mental assumption of bimodality, sample histogram fits 
for three representative magnitude ranges are shown in 
Figure [51 Two modes are definitely required, but the ar- 
bitrary addition of more than two subcomponents results 
in no for mal improvement to the total color distribution 
(see also lLarsen et al.|[2001t [Harris 2009; Cock croft et al.l 



120091 for similar tests in several other galaxies), showing 
that bimodality is both a necessary and sufficient condi- 
tion for describing the broadband-color distribution. 

7.1. The Blue Sequence 

The behavior of the key bimodal parameters (the 
means fii,fJ.2 and the dispersions u\,a2) with luminos- 
ity is shown in Figure [9l A nonzero MMR will appear as 
a systematic trend of the mean color /1 with magnitude. 
By using the bimodal fits to define the central colors of 
each mode, we also explicitly account for any overlap of 
the color distributions and possible asymmetric biasing 
due to contamination of one mode by the other. We 
concentrate first on the trends along the blue sequence, 
which to date is the only one of the two claimed to show 
a clear MMR, from the studies of many other galaxies. 

The blue sequence in M87 shows up as a narrow, rather 
sharply defined feature in the CMD that is easily traced 
over the entire spatial region studied here. The nar- 
row spread and large population on the blue sequence 
mean that the numerical solution for its position and 
shape is very stable. In the subsample used here for 
the fits, it takes up 55-60% of the total GC population, 
and the p\ ratio stays nearly constant at all magnitude 
levels (upper panel of Figure[l0|). 2 Comparison of the re- 
sults in Tables 2 and 3 (and displayed in Figure [9|) shows 
that no significant changes in [i\ or <j\ occur when the 
red sequence is arbitrarily required to fall at a constant 
color. The mean color (normalized to R = 1.0 R e ff) 
is (g — i)(blue) = 0.80 ± 0.01 over the fainter half of 
the i— magnitude range studied here. Its intrinsic width 
has a dispersion rri = 0.06 ± 0.01 that remains highly 
consistent with magnitude, except in the faintest range 
i > 22.5 where it is partially broadened by p hotomet- 
ric measurement scatter. The [Maraston (2005) popula- 
tion models used above then show that the dispersion in 
color corresponds to an intrinsic rms met allicity spread 
o-(Fe/H) = 0.30 ± 0.08. For comparison, iHarrid ([20091 ) 
found cr(Fe/H)(blue) = 0.27 for the average of six other 
BCGs, measured from the (B — I) color index. 

The bimodal fits show that the blue sequence be- 
comes smoothly and distinctly redder with increasing 
luminosity, especially for i < 21.5 (Mr < —10). Over 
this range, the mean color behaves as (g — i) nor rn = 
(0.845 ± 0.090) - (0.021 ± 0.004) (i - 20), a slope that 
is different from zero at the 5a level. The resultant scal- 
ing of GC heavy-element abundance wi th luminosity is 
then Z ~ L 25±0 - 08 . F or comparison, iCockcroft et al.1 
(2009) and|HarrisJ (|00l) find an average blue-sequence 
scaling of Z ~ £°- 3 ±o.05 f or more than a dozen large E 
galaxies measured from a variety of both ground-based 
and HST data and applying to a very similar GC lumi- 
nosity range, though these studies also indicate that real 
differences from one galaxy to the next may exist. 

At fainter levels (Mi > — 10, corresponding roughly to 
lO 6 Af ), it is not clear that any significant trend exists 
and that the blue sequence could better be considered 
vertical. In summary, these results for M87 are closely 

2 This pi value is only a local one for the spatial region used here. 
The global one over the entire M87 halo will be larger because the 
blue GCs extend much further out into the halo. For comparison, 
IPene et al.l (2008) find p 1 (global) = 0.73 for M87 by integrating 
the number counts for both types of clusters to large radius. 
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similar to those of lHarrid (|2009f) for six other BCG giants: 
the MMR for the blue sequence in M87 exists, but it is 
nonlinear in form and most strongly affects the most mas- 
sive clusters. Comparison of the Megacam result with the 
recent thor ough analys i s of t he deep (V — I) data from 
HST/ACS ([Peng et all 12009( 1 also shows close similari- 
ties: they found a clearly nonzero trend towards redder 
colors over the magnitude range Mj < —9 and a nearly 
vertical sequence below that. 

The data discussed here emphasize that detecting the 
MMR correlation on observational grounds depends not 
just on using a metallicity-sensitive photometric index, 
but also on having a very large sample that extends to 
high luminosity. If we had restricted our analysis to just 
the inner zone (see, for example, the first panel of Fig- 
ure [6]) , the MMR with its n onlinear form would be es- 
sentially undetectable. The iWaters et al.1 (|2009f l study, 
which uses (V — I) colors for the clusters that fall in a 
single HST/ACS field, covers only the region within 2' 
of the center of M87. Their conclusion that no MMR is 
present in fact does not disagree with the present Mega- 
cam results for two reasons: first, at lower GC luminosity 
Mi > — 10 where the vast majority of the clusters fall, 
the MMR actually does become more nearly vertical and 
no trend is expected. Second, the smaller sample size in 
the HST/ACS data means that only a few dozen of its 
clusters lie above the transition level Mj ~ —10 where 
the MMR starts to become significant, making it hard to 
distinguish above the intrinsic scatter in the sequence. 
A more detailed co mparative analysi s of the HST/ACS 
data is provided bv lPeng et alj (|2009l ). 

7.2. The Red Sequence 

The red sequence has a typical mean color (g-i)norm = 
1.07 averaged over all magnitudes, but it is less pre- 
cisely definable at any level because its intrinsic disper- 
sion <T2 — 0.10 ± 0.01 is nearly twice as broad as that 
of the blue sequence. This larger spread also makes it 
more vulnerable to the field contamination at the upper 
end. The corresponding metallicity spread is g(Fe/H) = 
0.50 ± 0.13. From their (C -T t ) photometry. iForte et all 
(|2007( 1 measure a ratio cri/cra = 0.60, in exact agreement 
with o u rs. Fo r a combination of six other BCG galaxies, 
IHarrisI (|20fH) found cr(Fe/H)(red) = 0.45 from (B - I), 
also in good agreement with the M87 sequences. 

In previo u s stu d ies with oth e r ga la xies (notably 
Harris et alj [2001 IMieske et all [2001 IHarrisI [20091 : 
Cockcroft et all l200a ~ lPeiig et all I2009L where explicit 
tests have been made) no significant MMR's have been 
claimed to exist along the red sequence. On theoretical 
grounds, if self-enrichment is the physical origin of the 
MMR then we do not expect to find any correlation of 
metallicity with mass along the red sequence except for 
extremely massive clusters (several million Solar masses 
and higher) where the amount of self-enrichment can fi- 
nall y reach above the alre ady-large pre-enrichment level 
(see lBailin fc Harrisll2009( ). In most galaxies neither the 
blue nor red sequences reach that high. However, the 
large and homogeneous database for M87 offers a chance 
to test this notion empirically. 

For the fainter range i > 21.5, no significant change 
shows up in the mean color n%. At brighter levels, the 
formal solutions show a surprisingly large change towards 
bluer mean color at higher luminosity; that is, the two se- 



quences give the appearance of converging towards each 
other there. Taken at face value, the red sequence would 
then display a "negative MMR" towards lower metallic- 
ity at higher mass, contrary to any existing model in- 
terpretation. Nominally, the slope for the unconstrained 
fits in Table 2 is A(g — i)/Ai ~ 0.03, corresponding to 
Z ~ L"(0-3±o.i5) 

Closer examination indicates, however, that this trend 
is not as robust as for the blue sequence. Formally the 
slope is significant at only the 2 — a level. Perhaps more 
importantly, the combined effects of its higher internal 
spread and the significant amount of field contamina- 
tion in the CMD make the mean position of the sparsely 
populated upper end of the red sequence quite uncer- 
tain. This is shown in Table 4 and the upper panel of 
Figure O The numerical tests described above show that 
constraining the position of the blue sequence /ii has a 
large effect on the deduced position of /12 , pulling it much 
further to the blue. By contrast, the opposite test of con- 
straining the red sequence color has no important effect 
on the blue sequence (lower pair of lines in Figure [9]). 

Another result of these tests is shown in 
the lower panel of Figure [TO] Here, the ratio 
X 2 (constrained)/x 2 (unconstrained) is plotted for 
the different sets of bimodal fits, to demonstrate what 
happens to the goodness-of-fit of each solution relative 
to the "baseline" model fit of Table 2. These % 2 ratios 
are listed in the last columns of Tables 3 and 4. If 
we arbitrarily fix the red sequence (Table 3, and red 
line in Figure [TO])), then the quality of fit remains 
virtually as good as the baseline model except for the 
very brightest bin. But for the opposite case where we 
fix the blue sequence (Table 4, and blue line in Figure 
HOP , a much more dramatic change occurs: at the bright 
end, the quality of fit degrades strongly. The same is 
true for the final set of fits where both sequences are 
arbitrarily constrained to their mean values. The case 
where the blue sequence is fixed also leads to a steep 
and unrealistic decrease in the p\ ratio at the bright end 
(blue line in the upper panel of Figure 1 10|) . 

In other words, over the range i > 21 [Mi > —10.5), 
the assumption that both sequences are strictly vertical 
(no MMR) provides an entirely reasonable fit to the data. 
Brighter than this, the null hypothesis of no MMR be- 
comes increasingly invalid, and closer attention must be 
paid to the curvature in one or both of the sequences. 

In summary, the present data point clearly to the exis- 
tence of a blue-sequence MMR. Correctly representing it 
is crucial for the quality of fit of the entire color distribu- 
tion. By contrast, the data indicate - though much more 
hesitantly - that a negative MMR may exist along the 
red sequence, though we cannot reject the null hypothe- 
sis (no MMR) as yet: the numerical tests described above 
show that ignoring the red MMR has very little effect on 
the quality of the bim odal fits. 

The recent study of iPeng et al . (2009) uses a sample 
of ~ 2000 M87 clusters measured in the less sensitive 
(V — I) color, but has the advantage that it is almost 
completely free of field contamination. They find that 
the red sequence has no significant slope at all luminosity 
levels. 

8. DISCUSSION AND INTERPRETATION 
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The results of this study fall into line with the most 
com plete recent co llections of data for o ther large galax- 
ies (jMieske et al.l 120061 lHarrisI 120091: iCockcroft et~aTl 
2009). Above a transition point of about 10 6 M Q , the 
low-metallicity clusters become progressively more en- 
riched at higher luminosity, scaling as Z ~ i 025 . 
This nonlinear MMR can be reasonabl y well understood 
through a basic self-enrichm ent model (jStrader fc Smith! 
120081: feailin fc Harrisl[2009l ). 

At the same time, a MMR of similar amplitude and 
reverse direction - but of much lower significance - is 
found alon g the red sequence, Z ~ £-(°- 3 ± - 15 ). in- 
terestingly, [Cockcroftirin (|2009f) find a mean scaling 
Z(red sequence) ~ L-i - 1 ^ - 1 ) for the average of 15 gi- 
ant E galaxies, a result that is only marginally different 
from zero but at least points in the same direction as 
the M87 data. Thus it may be worth considering how 
the red sequence could have embedded in it some form 
of substructure appearing as a MMR of its own in a way 
that does not require self-enrichment. 

To discuss the possibilities, we first ask what the es- 
sential differences are between the two modes. The 
blue, metal-poor clusters make up a narrow, well de- 
fined sequence that is replicated in very much the same 
way over and over in galaxies of all types, differing 
primarily just in their total numbers within different 
galaxies. Their mean metallicity (Fe/H) ~ —1.4 is 
only weakly correlated with galaxy mass and shows 
very much the same internal dis persion g(Fe/H) ~ 0.3 
from one system to another (e.g. IBurgarella et al.l 



Strader et aD 120041 : iBrodie fc Straderl 120061 : iPeng et all 
2006f ). In addition, the blue GCs display a consis- 



tently old ~ 13— Gy mean ag e with little internal age 
spread (e.g. iPuzia et al 1 120051 IBrodie fc Straderl 120061: 
IBeaslev et al.l 120081: iMarin-Franch et all l2009f ) . Many 
lines of argument advanced in recent years have built 
up the view that the metal-poor GCs are likely to have 
formed in the first brief round of star formation within 
the "pregalactic dwarfs" that began their star formation 
just at the beginning of the major stage of hierarchi- 
cal merging, and additional metal-poor GCs could have 
accreted l ater as part of their host dwarf satellites (see 
especially iSearle fc Zinnl 119781: IHarris fc Pudrita 119941: 
Forbes et al.l 119971: iCote et al.1 119981: IBurgarella et all 
20011: IBeaslev et al.l l2002t iKravtsov fc Gnedinl 120051 : 



Brodie fc Straderl 120061 : iBekki et alp 2008. among many 



others). These small host dwarfs with their shallow po- 
tential wells and initially metal-poor gas could support 
only one major round of star formation and thus left be- 
hind metal-poor stars plus much unconverted gas. 

By contrast, the red, metal-richer sequence is a 
more heterogeneous population that in some ways 
is harder to model in detail. Their mean metal- 
licity (Fe/H) ~ —0.4 correlates wea k ly with host 
galaxy luminosity dForbes et all 119971 : iStrader et alJ 
120041: IBrodie fc Straderl I2006D ~ Both their metallicity 
distribution and spatial distribution within the halos of 
their galaxies as a whole are similar to the field-star pop- 
ulation that has a broad, metal-rich MDF and makes up 
the bulk of a typical giant elliptical dGeisler et al.lfl~996t 
Harris fc Harris! [2002 iRhode fc Zep'J l2004t iForte et all 



chical mergi ng that built most of the ga laxy. Hierarchical 
models (e.g. IBeaslev et al.ll2002t 120031) show that multi- 
ple "in situ" starbursts, major mergers, and accretions 
of gas-rich satellites would all have contributed to this 
stage. If this basic picture is correct, then the final pop- 
ulation of metal-rich GCs would be expected to have a 
wide range of ages and metallicities, and more so in big- 
ger galaxies that experienced more starbursts and merg- 
ing events. Furthermore, the red-GC populations should 
also differ in detail from one galaxy to another depending 
on their individual merger and starburst histories. 

On the observational side, the evidence points in- 
creasingly to the same view that the metal-rich GC 
"population" is actually a composite of objects with 
various ages, with its subcomponent MDFs heavily over- 
lapped and thus of nearly similar mean metallicities (e.g . 
Puzia et al.ll2005l:lKundu et alJl2005b IBrodie et al.ll2005t 
Dirsch et al. 2005; Larsen et al. 2005; Brodie & Strader 



20061: ICenarro et al.l 120071: IGoudfrooii et all I2007L 
Beaslev et al.1 120081: iMavbhate et al.1 120091 ). Consistent 



2007). That is, the metal-rich GCs can reasonably be 



thought of as forming during the major stage of hierar- 



with this view, both the relative numbers and luminosity 
functions of the red sequence differ noticeably from one 
galaxy to another, so that the high-luminosity end of 
the red GCLF does not extend to the same maximum 
even in galaxies with similar total GC populations. For 
example, in NGC 1399 the blue and red sequences have 
very similar "upper ends" while in M87 the blue se- 
quence extends much higher (jForte et al. 2007]); see also 
the di fference between th e Antlia giants NGC 3258 and 
3268 (B assino et al.ll2008l ). By contrast with most other 
systems, in the Hydra and Coma supergiants NGC 3311 
and NGC 4874 the red sequence is the one that extends 
to higher luminosity, up to cluster mass e s ~ 10 7 M^ 
and perhaps higher (jWehner et alJ 120081 : IHarris et alJ 
2009). Observations of young, ma ssive star clusters 
in present-day st arburst galaxies (jLarsen fc Richtlerl 
2000; Larsen 2009) also show that the cluster formation 
efficiency and the initial cluster mass function can vary 
with the total star formation rate. 

These various pieces of evidence point to a way to un- 
derstand, at least in principle, how we might obtain a 
red GC sequence with an observable mass/metallicity re- 
lation that could be either positive or negative. In the 
simplest such case, suppose that the metal-rich GCs com- 
prise two main subpopulations that can be thought of as 
the products of the two biggest mergers or starbursts 
that assembled the main body of the galaxy. In addi- 
tion, suppose that the earlier and slightly less metal-rich 
of these two starbursts generated its GCs with a broader 
mass function that extended to higher luminosity than 
did the later and metal-richer burst. What we would 
see today in the GCS as a whole would then be a red se- 
quence with a broad MDF and a "negative MMR" at the 
top end. This population-based approach has nothing 
to do with any self-enrichment, is completely consistent 
with the evidence summarized above, and does not in- 
voke any unusual physical effect. Simulated GCSs along 
these lines are now being constructed to test these ideas 
further and will be discussed in a later paper. 

9. SUMMARY 

Wide-field photometry with CFHT/Megacam in gri 
has been used to study the properties of the M87 globu- 
lar cluster system. Approximately 6200 GCs have been 
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measured to a limit i ~ 23.0 (equivalent roughly to 
Mj ~ —8.5, a little less than half the total population). 
The primary results of this study are these: 

1. The metal-poor (blue) and metal-rich (red) clus- 
ter sequences are clearly delineated especially in 
the (g — i) color index. While the red sequence 
does not extend detectably outward to projected 
galactocentric distances beyond ^50 kpc against 
the field contamination, the blue sequence can be 
seen to 100 kpc and beyond. 

2. Both the red and blue subpopulations have shallow 
but significant radial metallicity gradients, corre- 
sponding to heavy-element abundance scaling with 
distance of Z ~ R~°- 12 (blue) and Z ~ R-° 17 
(red). 

3. The blue sequence shows a nonlinear 
mass/metallicity relation: GC metallicity in- 
creases with luminosity. This relation is nearly 
vertical (that is, mean color is uncorrelated with 
metallicity) for Mr > — 10, while at brighter levels 
we find a trend roughly reproduced by Z ~ L 25 . 
This result agrees well with the findings in nu- 
merous other giant E galaxies, as well as the best 
current data on M87 itself. Finding the trend 
depends critically on being able to sample the 
GC population to its highest luminosity levels. 
The blue-sequence MMR can be physically under- 
stood as the result of GC self-enrichment during 
formation. 



4. The red sequence is almost twice as broad as the 
blue sequence in its internal scatter of colors. The 
formal numerical solutions show that it has a "neg- 
ative MMR" at the top end that is at least as 
steep as the blue-sequence MMR. However, the 
same numerical solutions show that this trend is 
not strongly significant, and simply assuming a 
constant mean color for the red sequence does not 
noticeably decrease the quality of the bimodal fits. 

5. In giant galaxies like M87, the red-GC sequence 
is likely to be a composite population, unlike the 
more internally homogeneous blue sequence. Every 
major starburst or gas- rich merger that contributed 
to the formation of the main body of the galaxy 
should have contributed to its metal-rich GC pop- 
ulation. If the individual major starbursts gener- 
ated GCs following slightly different initial mass 
functions, then features such as the negative-MMR 
can be built in to the total red sequence that we 
see today. 
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TABLE 1 
CFHT/Megacam Images 



Filter JV(cxp) Total t (sec) FWHM 



9' 


11 


3520 


0.74" 


r' 


10 


3000 


0.80" 


i' 


11 


3300 


0.82" 



TABLE 2 
Radial Distributions 



R Range R n c i a c i (arcmin 2 ) 



Blue GCs 








1.00' 


- 1.30' 


1.14' 


68.4 


31.57 ±3.83 


1.30 


- 1.69 


1.48 


84.1 


22.94 ±2.52 


1.69 


- 2.20 


1.93 


128.4 


20.75 ± 1.84 


2.20 


- 2.86 


2.50 


135.3 


12.93 ± 1.13 


2.86 


- 3.71 


3.26 


206.4 


11.67 ±0.83 


3.71 


-4.83 


4.23 


256.2 


8.57 ±0.55 


4.83 


- 6.27 


5.50 


293.9 


5.82 ±0.36 


6.27 


- 8.16 


7.15 


388.8 


4.56 ±0.26 


8.16 - 


- 10.60 


9.30 


428.5 


2.97 ±0.18 


10.60 


- 13.79 


12.09 


423.6 


1.74 ±0.14 


13.79 


- 17.92 


15.72 


426.8 


1.04 ±0.12 


17.92 


- 23.30 


20.43 


472.9 


0.68 ±0.11 


23.30 


- 30.29 


26.56 


382.1 


0.33 ±0.10 


Red GCs 








1.00' 


- 1.30' 


1.14' 


81.2 


37.47 ±4.18 


1.30 


- 1.69 


1.48 


102.7 


28.03 ± 2.79 


1.69 


- 2.20 


1.93 


135.8 


21.93 ± 1.90 


2.20 


- 2.86 


2.50 


143.2 


13.69 ± 1.16 


2.86 


- 3.71 


3.26 


167.6 


9.48 ± 0.75 


3.71 


-4.83 


4.23 


170.2 


5.70 ±0.46 


4.83 


-6.27 


5.50 


186.8 


3.70 ± 0.30 


6.27 


- 8.16 


7.15 


189.3 


2.22 ±0.20 


8.16 - 


- 10.60 


9.30 


187.1 


1.30 ±0.15 


10.60 


- 13.79 


12.09 


166.2 


0.68 ±0.12 


13.79 


- 17.92 


15.72 


145.7 


0.35 ±0.11 


17.92 


- 23.30 


20.43 


94.4 


0.14 ±0.10 
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TABLE 3 

RMIX Bimodal Fits to the (g - i) Colors 



i range JV(bin) Ml (±) <*\ (±) (±) &2 (±) Pi (±) 



(23.00,23.25) 
(22.75, 23.00) 
(22.50,22.75) 
(22.25,22.50) 
(22.00,22.25) 
(21.75,22.00) 
(21.50, 21.75) 
(21.25,21.50) 
(21.00,21.25) 
(20.75,21.00) 
(20.50, 20.75) 
(20.00, 20.50) 
(19.50, 20.00) 



455 


0. 


497 


0. 


486 


0. 


397 


0. 


330 


0. 


311 


0. 


257 


0. 


191 


0. 


152 


0. 


138 


0. 


115 


0. 


122 


0. 


57 


0. 



809(0.009) 
798(0.008) 
799(0.006) 
790(0.007) 
807(0.007) 
802(0.006) 
806(0.008) 
818(0.012) 
829(0.010) 
838(0.009) 
823(0.009) 
846(0.010) 
844(0.017) 



0.085(0.007) 
0.079(0.006) 
0.064(0.004) 
0.066(0.006) 
0.064(0.005) 
0.059(0.005) 
0.061(0.006) 
0.068(0.008) 
0.052(0.011) 
0.060(0.000) 
0.060(0.000) 
0.060(0.000) 
0.060(0.000) 



1.097(0.018) 
1.097(0.013) 
1.078(0.013) 
1.079(0.017) 
1.083(0.013) 
1.077(0.013) 
1.073(0.022) 
1.067(0.030) 
0.996(0.040) 
1.051(0.020) 
1.026(0.021) 
1.033(0.021) 
0.985(0.033) 



0.099(0.013) 
0.092(0.010) 
0.100(0.009) 
0.116(0.013) 
0.086(0.010) 
0.098(0.010) 
0.105(0.015) 
0.095(0.019) 
0.122(0.019) 
0.100(0.000) 
0.100(0.000) 
0.100(0.000) 
0.100(0.000) 



0.633(0.035) 
0.616(0.035) 
0.576(0.035) 
0.542(0.035) 
0.606(0.035) 
0.553(0.035) 
0.592(0.035) 
0.625(0.035) 
0.415(0.035) 
0.591(0.058) 
0.586(0.064) 
0.540(0.072) 
0.466(0.143) 



TABLE 4 

Bimodal Fits (Red Sequence Constrained) 



i range JV(bin) /hi (±) a\ (±) <T2 (±) pi (±) \l(con) /xl(uncon) 



(23.00, 
(22.75, 
(22.50, 
(22.25, 
(22.00, 
(21.75, 
(21.50, 
(21.25, 
(21.00, 
(20.75, 
(20.50, 
(20.00, 
(19.50, 



23.25) 
23.00) 
22.75) 
22.50) 
22.25) 
22.00) 
21.75) 
21.50) 
21.25) 
21.00) 
20.75) 
20.50) 
20.00) 



455 
497 
486 
397 
330 
311 
257 
191 
152 
138 
115 
122 
57 



0.799(0.009) 
0.788(0.006) 
0.797(0.005) 
0.787(0.006) 
0.803(0.006) 
0.800(0.005) 
0.805(0.006) 
0.819(0.008) 
0.848(0.010) 
0.841(0.008) 
0.829(0.008) 
0.853(0.009) 
0.871(0.011) 



0.079(0.005) 
0.072(0.005) 
0.062(0.004) 
0.064(0.005) 
0.061(0.004) 
0.057(0.004) 
0.060(0.005) 
0.068(0.006) 
0.068(0.008) 
0.060(0.000) 
0.060(0.000) 
0.060(0.000) 
0.060(0.000) 



0.114(0.009) 
0.108(0.007) 
0.104(0.007) 
0.121(0.009) 
0.093(0.008) 
0.102(0.008) 
0.106(0.010) 
0.093(0.011) 
0.094(0.013) 
0.100(0.000) 
0.100(0.000) 
0.100(0.000) 
0.100(0.000) 



0.578(0.030) 
0.562(0.028) 
0.562(0.027) 
0.527(0.031) 
0.582(0.032) 
0.542(0.033) 
0.587(0.038) 
0.631(0.043) 
0.676(0.056) 
0.618(0.048) 
0.640(0.051) 
0.605(0.054) 
0.775(0.074) 



1.08 
1.11 
0.90 
0.92 
0.91 
0.92 
0.88 
0.88 
1.11 
0.99 
1.13 
1.25 
1.82 



TABLE 5 

Bimodal Fits (Blue Sequence Constrained) 



i range N(bin) a\ (±) ^ 2 (±) &2 (±) Pi (±) xl(con)/xl(uncon) 



(23.00, 
(22.75, 
(22.50, 
(22.25, 
(22.00, 
(21.75, 
(21.50, 
(21.25, 
(21.00, 
(20.75, 
(20.50, 
(20.00, 
(19.50, 



23.25) 
23.00) 
22.75) 
22.50) 
22.25) 
22.00) 
21.75) 
21.50) 
21.25) 
21.00) 
20.75) 
20.50) 
20.00) 



455 
497 
486 
397 
330 
311 
257 
191 
152 
138 
115 
122 
57 



0.081(0.005) 
0.080(0.005) 
0.064(0.004) 
0.071(0.005) 
0.061(0.005) 
0.058(0.004) 
0.058(0.005) 
0.059(0.007) 
0.042(0.009) 
0.060(0.000) 
0.060(0.000) 
0.060(0.000) 
0.060(0.000) 



1.083(0.014) 
1.100(0.010) 
1.079(0.010) 
1.092(0.012) 
1.074(0.012) 
1.076(0.012) 
1.063(0.018) 
1.023(0.027) 
0.965(0.019) 
1.011(0.017) 
1.004(0.020) 
0.990(0.018) 
0.918(0.028) 



0.108(0.011) 
0.090(0.008) 
0.099(0.008) 
0.107(0.009) 
0.093(0.009) 
0.099(0.009) 
0.111(0.013) 
0.119(0.017) 
0.123(0.010) 
0.100(0.000) 

o.ioo(o.ooo) 

0.100(0.000) 

o.ioo(o.ooo) 



0.600(0.034) 
0.621(0.027) 
0.579(0.028) 
0.575(0.031) 
0.582(0.035) 
0.549(0.034) 
0.569(0.044) 
0.506(0.074) 
0.301(0.083) 
0.451(0.061) 
0.513(0.066) 
0.347(0.076) 
0.063(0.206) 



0.98 
0.88 
0.88 
1.22 
0.92 
0.88 
0.91 
0.98 
1.36 
1.49 
1.16 
6.83 
1.83 
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Fig. 1. — Use of SExtractor parameters for rejection of nonstellar objects. The measured FWHM and effective radius of each detected 
object on the M87 field are plotted against aperture magnitude. Objects lying outside either of the marked boxes are rejected from 
the detection lists (see text). The example shown here is from the r— band frame; note that the r values plotted here are only the raw 
instrumental magnitudes and should not be compared with the final color-magnitude diagrams below. 
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Fig. 2. — Residuals in the measured positions for the 932 stars in the USNO-B1 catalog that were found to lie within the M87 Mcgacam 
field, plotted as a function of declination in degrees. The rms scatter around zero is ±0.23" in both right ascension and declination. 
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Fig. 3. — Final color-magnitude diagrams for the starlike and near-starlike objects detected in all three frames (g,r,i). Objects lying 
outside either of the marked boxes are rejected as field contamination (mostly foreground Milky Way halo or disk stars). 



16 Harris 







10 



R (arcmin) 
20 30 



40 







I — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | 
50 100 150 kpc 



i 1 r 



~i 1 1 1 1 r 



~i 1 1 1 r 



1.4 



1.2 4: 



0.8 



0.6 



j i i i i i i i i i i i i i i i i i i i i i i_ 







10 15 
R/R eff 



20 



25 



Fig. 4. — Measured (g — i) color index versus projected galactocentric distance for the "starlike" objects brighter than i = 23.0. This 
sample is dominated by globular clusters (see text). Radii arc plotted in units of the M87 spheroid effective radius (bottom axis) or in 
kiloparsecs or arcminutes (top axes). The calculated radial metallicity gradients for the blue and red GC sequences are shown by the solid 
lines. 
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Fig. 5. — Spatial distributions of GC candidate objects brighter than i = 23.0 as in the previous figure, subdivided into the blue sequence 
(left panel) and red sequence (right panel). The much larger spatial extent of the blue GC population is evident. Note also the two vacant 
horizontal stripes above and below center that represent the large gaps between the ranks of CCDs in the camera mosaic. 



18 



Harris 



19 - 



20 



21 



22 



23 



24 - 



1'<R<3' 



iii|iiiiiiiiiii 



-SAHv v.- 



i i i i i i i i 



0.8 1.2 
g-i 



3'-6' 



iii|iiiiiiiiiii 



6'-15' 



15'-27' 




1 1 1 1 1 1 



0.8 1.2 
g-i 




27' < R 



i.i i ti 1 1 1 1 1 1 1 1 1 1 



W 'Crt .:•: ••»:'• - 
.... 




: "/ — 



I I I I I I I I 



0.8 1.2 
g-i 



Fig. 6. — Color-magnitude data in (i, g — i) for the starlike objects, separated by radial zone in projected galactocentric distance R. Five 
zones are shown. 
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Fig. 7. — Projected radial profiles for the blue and red GC subsystems, using the i < 23.0 data from Figure [4] The number density of 
GCs, a c i = <t — crbkgd, is plotted in power-law form (left panel) and again in the de Vaucouleurs-law form (right panel). The blue-sequence 
GCs are shown as solid dots and the red-sequence GCs as open triangles. In the left panel the Hubble profiles fitting each of the two 
components as given in the text are superimposed. 
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Fig. 8. — Sample RMIX solutions for the (g — i) color distributions in three one-magnitude bins in i as labelled. The (g — i) colors of the 
individual clusters are normalized to a galactocentric distance R = l.0R e ff to remove the spatial metallicity gradients in the system (see 
text). In each panel, the dashed lines show the Gaussian curves matching the blue and red sequences, while the solid lines show the sum 
of the two components. 
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Fig. 9. — Upper panel: Binned mean colors in (g — i) for the blue (solid dots) and red (open triangles) sequences, as listed in Table 
\3\ Lower panel: Internal standard deviation (color width) of the blue and red sequences as a function of i magnitude. For i < 21, the 
a— values are constrained at their mean levels (see text). For the blue sequence (lower curve in each panel), the dashed lines show how the 
fitted color and dispersion change if the red sequence is constrained to have fi2 = 1.07 at all magnitudes, as listed in Table 3. For the red 
sequence (upper curve in each panel), the dashed lines show the changes if the blue sequence is constrained to have color p^i = 0.80 at all 
magnitudes. 
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Fig. 10. — Upper panel: The fraction pi of the total GC population in the blue sequence, plotted as a function of i magnitude and obtained 
from the bimodal RMIX fits. The solid line shows the results from Table 2, where both blue and red means p-1,1^2 are unconstrained. The 
red long-dashed line shows the results where the red sequence is constrained to lie at pi2 = 1-070; the blue short-dashed line shows the 
results where the blue sequence is constrained at /ii = 0.80; and the dotted line shows the results where both sequences are constrained. 
Lower panel: The quality-of-fit ratio \ t (constrained) /x^ (unconstrained) for three different versions of the bimodal-sequence fits. The solid 
line shows the ratio for the case where both sequences are constrained at fii = 0.80, fi2 = 1.07; the red long-dashed line shows the case for 
p.2 = 1.07; and the blue short-dashed line shows the case for = 0.80. A ratio larger than 1.00 means that the solution provides a worse 
fit than the unconstrained case. 
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Fig. 11. — Color-magnitude diagram in [i, g — i) for the data used in the bimodal sequence fits (see text). The solid lines connect the 
mean points from Table [3] from the unconstrained bimodal fits. 



